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Abstract 



We calculate the critical value of the hopping parameter, Kc, in Lattice 
QCD with Wilson fermions, to two loops in perturbation theory. 

This quantity is an additive renormalization; as such, it is characterized 
not only by the standard caveats regarding the asymptotic nature of pertur- 
bative results, but also by a linear divergence in the lattice spacing. Conse- 
quently, our calculation tests rather stringently the limits of applicability of 
perturbation theory. 

We compare our results to non perturbative evaluations of Kc coming from 
Monte Carlo simulations. 

Finally, we apply a tadpole improvement technique on our results; this 
shifts them quite favourably towards the non perturbative values. 
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I. INTRODUCTION 



In this paper we study the hopping parameter in lattice QCD with Wilson fermions. In 
particular, we compute its critical value to two loops in perturbation theory. 

Wilson fermions are the most straightforward and widely used implementation of 
fermionic actions on the lattice. This implementation circumvents the fermion doubling 
problem by introducing a higher derivative term with a vanishing classical continuum limit, 
to lift unphysical propagator poles completely. At the same time, the action is strictly local, 
which is very advantageous for numerical simulation. 

The price one pays for strict locality and absence of doublers is, of course, well known: 
The higher derivative term breaks chiral invariance explicitly. Thus, merely setting the 
bare fermionic mass to zero is not sufficient to ensure chiral symmetry in the quantum con- 
tinuum limit; quantum corrections introduce an additive renormalization to the fermionic 
mass, which must then be fine tuned to have a vanishing renormalized value. Consequently, 
the hopping parameter k, which is very simply related to the fermion mass, must be appro- 
priately shifted from its naive value, to recover chiral invariance. 

By dimensional power counting, the additive mass renormalization is seen to be linearly 
divergent with the lattice spacing. This adverse feature of Wilson fermions poses an addi- 
tional problem to a perturbative treatment, aside from the usual issues related to lack of 
Borel summability. Indeed, our calculation serves as a check on the limits of applicability of 
perturbation theory, by comparison with non perturbative results coming from Monte Carlo 
simulations. 

Starting from our two-loop results, we also provide improved estimates of the critical 
value of by performing a resummation to all orders of cactus diagrams These are 
tadpole-like diagrams which are gauge invariant and dress the propagators and vertices in 
our calculation. This improvement technique, among others, has so far been applied mostly 
to the one-loop multiplicative renormalization of various operators . It is interesting to 
explore to what extent such methods lead to an improvement even in a sensitive case such 
as the one at hand. We find that our improved estimates compare quite well with Monte 
Carlo data also in this case. 

The paper is organized as follows: In Sec. |T| we define the quantities which we set out to 
compute, and describe our calculation. In Sec. |IT| we present our results and compare with 
Monte Carlo evaluations. In Sec. |^ we obtain the improved estimates coming from cactus 
resummation. 



II. FORMULATION OF THE PROBLEM 

QCD with Wilson fermions on the lattice is described by the following action (see, e.g., 
Ref. m for standard notation and conventions): 

= - E Tr [1 - U,,{x)]+Y.Y.U^)D{x,y)Uy) (1) 

fi'O x,fj,,u i=l x,y 

U^u{x) is the standard product of link variables U^^y around a plaquette in the direction 
/i— z/, originating at point x, and D{x,y) is given by: 
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D{X, y) — aniB S^^y + ^ X! t^M {Ux,ySx+fl,y — Ux,y5x,y+[l) — r {Ux,y5x+[l,y — '2Sx,y + Ux,ySx,y+fi)] 

(2) 

As usual, go denotes the bare coupling constant and a is the lattice spacing. The bare 
fermionic mass uib must be set to zero for chiral invariance in the classical continuum limit. 

The higher derivative term, multiplied by the Wilson coefficient r, breaks chiral in- 
variance. It vanishes in the classical continuum limit; at the quantum level, it induces 
nonvanishing, flavor-independent corrections to the fermion masses. 

Numerical simulation algorithms usually employ the hopping parameter, 

K = (3) 

as a tunable quantity. Its critical value, at which chiral symmetry is restored, is thus l/8r 
classically, but gets shifted by quantum effects. 

The renormalized mass can be calculated in textbook fashion from the fermion self- 
energy. Denoting by 'E^{p,mB, go) the truncated, one particle irreducible fermionic two- 
point function, we have for the fermionic propagator: 

1 °° f 1 \'' 

'^(^^) ^ ^Y~Z — TtI^ ^^{P:'mB,go)TY- — rr 
If + m{p) ^ \ If + m{p) J 

= [i^ + m{p) - E^{p, TUB, go) (4) 
where : = ^ 7^- sin(ap'*), m{p) — uib H y^sin^(ap^/2). 

Requiring that the renormalized mass vanish, leads to: 

5-1(0) = =^ ms = E^(0,mB,^o) (5) 

The above is a recursive equation for ms, which can be solved order by order in perturbation 
theory. 

We write the loop expansion of as: 

E^(0,mB,^o)=^o'S^'^+^oS('^ + --- (6) 

Fig. I shows the two diagrams contributing to the 1-loop result E*^i\ The fermion mass 
involved in these diagrams must be set to its tree level value, uib — > 0. The i*^ diagram 
gives a contribution of the form cf'\ where 62'^ are numerical constants. 
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FIGURE I. One-loop diagrams contributing to S^. Wavy (solid) lines represent gluons (fermions). 



A total of 26 diagrams contribute to the 2-loop quantity E^^^, shown in Fig. II. Genuine 
2-loop diagrams must again be evaluated at — >■ 0; in addition, one must include to 
this order the 1-loop diagram containing an 0{gQ) mass counterterm (diagram 23). The 
contribution of each diagram can be written in the form 

(^^-i)(4!' + f + f4?) (7) 

where c^j , , cfj are numerical constants. Certain sets of diagrams, corresponding to renor- 
malization of loop propagators, must be evaluated together in order to obtain an infrared- 
convergent result: these are diagrams 7+8+9+10+11, 12+13, 14+15+16+17+18, 19+20, 
21+22+23. 
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FIGURE II. Two- loop diagrams contributing to S^. Wavy (solid, dotted) lines represent gluons 
(fermions, ghosts). Crosses denote vertices stemming from the measure part of the action; a solid 
circle is a fermion mass counterterm. 
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III. NUMERICAL RESULTS 



The evaluation of the diagrams in this computation requires very extensive analytical 
work. To this end, we use a Mathematica package which we have developed for symbolic 
manipulations in lattice perturbation theory (see, e.g., Ref. 0). Applied to the present 
case, this package allows us to perform in a rather straightforward way the following tasks: 
Contraction among the appropriate vertices; simplification of color/Dirac matrices; use of 
trigonometry and momentum symmetries for reduction to a more compact, canonical form; 
automatic generation of highly optimized Fortran code for the loop integration of each type 
of expression. 

The integrals, typically consisting of a sum over a few hundred trigonometric products, 
are then performed numerically on lattices of varying finite size L. Our programs per- 
form extrapolations of each expression to a broad spectrum of functional forms of the type: 
ejj(lnL)-^/L*, analyze the quality of each extrapolation using a variety of criteria and 
assign statistical weights to them, and finally produce a quite reliable estimate of the sys- 
tematic error. Taking L < 28 leads to a sufficient number of significant digits in the results 
we present. 

One important consistency check can be performed on those diagrams which are sepa- 
rately IR divergent; taken together in groups, as listed below Eq. (|^), they give a finite and 
very stable extrapolation. 

fl) (2) (2) (2) 

We present below the numerical values of the constants q , ^ , c^^, c\ i. These constants 
depend only on the Wilson parameter r; following common practice, we set r = 1. 

Table | contains the contributions to the 1-loop quantity The total 1-loop result is 

iV^ - 1 

= ^ (-0.162857058711(2)) (8) 

This result is known in the literature (see, e.g., Ref. 0, p. 246, and references contained 
therein). 

The contributions to the 2-loop quantity S*^^^ are presented in table |T[ The total 2-loop 
result is 

= (Ar2 _ 1)(-0.017537(3) + ^ 0.016567(2) + ^ 0.00118618(8)) (9) 

In order to make a comparison with numerical simulations, let us set = 3, Nf = 2 in 
the above; we obtain 

(iV=3, Nf=2) = -0.434285489897(5) (10) 



E(2)(Ar=3,A^j=2) = -0.11925(3) (11) 

In Table |T| we compare the final results for m^^^ = g^T.^^^ and m^^^ = (^o^*-^-* + fi'o^'-^'' 
with numerical simulation data at values of /5 = G/g^ equal to 5.6 and 5.5 [§]. Also 
included in this table are the improved results obtained with the method described in the 
following Section. For easier reference. Table |I^ presents our results in terms of the critical 
hopping parameter k,c = 1/(2 + 8r). 
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IV. IMPROVED PERTURBATION THEORY 



In order to obtain improved estimates from lattice perturbation theory, one may perform 
a resummation to all orders of the so-called "cactus" diagrams jlj-Q . Briefly stated, these are 
gauge-invariant tadpole diagrams which become disconnected if any one of their vertices is 
removed. The original motivation of this procedure is the well known observation of "tadpole 
dominance" in lattice perturbation theory. In the following we refer to Ref. |l| for definitions 
and analytical results. 

Since the contribution of standard tadpole diagrams is not gauge invariant, the class 
of gauge invariant diagrams we are considering needs further specification. By the Baker- 
Campbell-Hausdorff (BCH) formula, the product of link variables along the perimeter of a 
plaquette can be written as 

jj — gigoAx.^i ^igoAx+fi,,^ ^-igoAx+,^,^t ^-igoAx^^ 

= exp [igo{A^,f, + A.^+^^^ - v4^+^,^ - A^^^) + O(fi'o)} 

= exp {ig,F^l + iglF^, + Oig^) } (12) 

The diagrams that we propose to resum to all orders are the cactus diagrams made of 
vertices containing -F^^^j, . Terms of this type come from the pure gluon part of the lattice 
action. These diagrams dress the transverse gluon propagator Pa leading to an improved 
propagator Pa \ which is a multiple of the bare transverse one: 

1 - w(^o) 

where the factor w{gQ) will depend on g^ and A^, but not on the momentum. The function 
w{go) can be extracted by an appropriate algebraic equation that has been derived in Ref. |I| 
and that can be easily solved numerically; for SU{3), w{go) satisfies: 



nV3 - 4n + 81 = 2gl u{g^) = ^\ (14) 



The vertices coming from the gluon part of the action, Eq. (|1|), get also dressed using a 
procedure similar to the one leading to Eq. ([T^) [|l|. Vertices coming from the fermionic 
action stay unchanged, since their definition contains no plaquettes on which to apply the 
linear BCH formula. 

One can apply the resummation of cactus diagrams to the calculation of additive and 
multiplicative renormalizations of lattice operators. Applied to a number of cases of in- 
terest this procedure yields remarkable improvements when compared with the avail- 
able nonperturbative estimates. As regards numerical comparison with other improvement 
schemes, such as boosted perturbation theory P,p!0|, cactus resummation fares equally well 
on all the cases studied ||^. 

One advantageous feature of cactus resummation, in comparison to other schemes of 
improved perturbation theory, is the possibility of systematically incorporating higher loop 
diagrams. The present calculation best exemplifies this feature, as we will now show. 
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Dressing the 1-loop result is quite straightforward: the fermionic propagator and vertices 
stay unchanged, and only the gluon propagator gets simply multiplied by 1/(1 — w{go). The 
resulting values: ^^c^dressed ^c^dressed; showu in Tables m and respectively. It 
is worth noting that these values already fare better than the much more cumbersome 
undressed 2- loop results. 

We now turn to dressing the 2-loop results. Here, one must take care to avoid double 
counting: A part of diagrams 7 and 14 has already been included in dressing the 1-loop 
result, and must be explicitly subtracted from S*^^-* before dressing. Fortunately, this part 
(we shall denote it by Sg^j^) is easy to identify, as it necessarily includes all of the 1/iV^ part 
m S(2). A simple exercise in contraction of SU (N) generators shows that E^Jj^ is proportional 
to {2N^ - 3){N^ - 1)/{3N^). There follows immediately that: 

= -0.016567(2A^2 _ ^-^^^^2 _ i)/(3Ar2) (15) 

(cf. Eq. I). 

A further complication is presented by gluon vertices. While the 3-gluon vertex dresses 
by a mere factor of (1 — w{gQ)), the dressed 4-gluon vertex contains a term which is not 
simply a multiple of its bare counterpart (see Appendix C of Ref. [0]). Once again, however, 
we are fortunate: this term must be dropped, being precisely the one which has already 
been taken into account in dressing the 1-loop result, while the remainder dresses in the 
same way as the 3-gluon vertex. In conclusion, cactus resummation applied to the 2-loop 
quantity E^^^ leads to the following rather simple recipe: 

™(2) _ ^(1) 90 ^y{2) _ y{2) ^ 9o 

'"c, dressed ~ ^ i /' ^ ^ V ^subJ ri / xio K^'^) 

,(2) _J ,.(2) 



The end results, dressed k), dressed 5 3.re included in Tables |T| and Comparing 
with the Monte Carlo estimates, we see a definite improvement over non-dressed values. At 
the same time, a sizeable discrepancy still remains, as was expected from start. This dis- 
crepancy sets a benchmark for lattice perturbation theory; multiplicative renormalizations, 
calculated to the same order and improved by cactus dressing, are expected to be much 
closer to their exact values. We hope to return to these calculations in a future publication. 
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TABLE 1. Coefficients cfK r = 1. 

i 4" 

1 -0.15493339023106 

2 -0.007923668480(2) 



TABLE n. Coefficients cf], cf]. r = l. 



.(2) 



-2.8 



.(2) 



3 0.002000362950707492 

4 0.00040921361(1) 

5 

6 -0.0000488891(8) 
7+8+9+10+11 -0.013927(3) 

12+13 
14+15+16+17+18 -0.005753(1) 

19+20 

21+22+23 0.000096768(4) 

24 

25 0.00007762(1) 

26 -0.00040000(5) 

27 

28 0.0000078-182(5) 



-0.00030005444260612375 

-0.00061382041(2) 



0.000097778(2) 

0.014525(2) 

0.00079263(8) 

0.0058323(7) 

0.000393556(7) 

-0.000096768(4) 



-0.00015524(3) 



-0.000006522(1) 

-0.000015090(1) 



TABLE IIL mi^^ and mi\ 



N = 3,Nf = 2, 



1. 





mi; 

(2) 

rric 
(1) 

c, dressed 

(2) 

c, dressed 

Simulation 



(3 = 5.5 



-0.473765988978(5) 

-0.61568(3) 

-0.658392964276(7) 

-0.76323(6) 

-0.8975 



(3 = 5.6 



-0.465305882032(5) 

-0.60219(3) 

-0.640695803036(7) 

-0.73997(6) 

-0.8446 



TABLE IV. k'^^ and k'^"' . N = 3, Nf = 2, r = 1. 



(2) AT- 



rvc 
^(2) 



K 



(1) 

c, dressed 



c, dressed 

Simulation 



J = 



0.1417943331149(4) 

0.147740(3) 

0.1496286052353(6) 

0.154475(6) 
0.16116(15) 



J = 5.6 



0.1414549557367(4) 

0.147154(3) 

0.1488403462990(6) 

0.153373(6) 
0.1585071^^ 
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